Automatic physical phasing X-ray crystallography 
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Phase invariants are important pieces of information about the atomic structures of crystals. 
There are several mathematical methods in X-ray crystallography to estimate phase invariants. 
The multi-wave diffraction phenomenon offers a unique opportunity of physically measuring phase 
invariants. In this work, the underneath principals for developing an automatic procedure to extract 
accurate phase-invariant values are described. A general systematic procedure is demonstrated, in 
practice, by analyzing intensity data from a KDP crystal. 
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I. INTRODUCTION 

In X-ray crystallography, the phases of the diffracted 
waves are roughly estimated by mathematical methods, 
know as Direct Methods^''^ , for analyzing intensity data 
sets composed of a large number of reflections. These 
methods exploit algebraic or probabilistic relationships 
among the phase values. Some of such relationships 
are triplet phase invariants; they are invariant from the 
choice of origin in the crystal lattice. Experimental pro- 
cedures allowing physical measurements of phase invari- 
ants are of great interest since, in principle, they could 
extend the efficiency of the Direct Methods to complex 
structures such as proteins. It would have to be com- 
pared to other procedures that are actually used to the 
same purposes, such as multiple anomalous dispersion 
and multiple isomorphous replacement'^. 

Physical measurements of triplet phase invariants are 
possible by means of three-beam diffraction (3BD) ex- 
perimental where the interference of simultaneously 
diffracted waves provide information on phase values. 
However, besides all experimental and analytical difh- 
culties involved in phase determination from 3BD ex- 
periments, the most serious and practical limitation of 
physical-phasing crystallography (PPC) is the reduce 
number of 3BD cases suitable for phasing. The reciprocal 
space of complex molecule crystals are full of reflections 
where isolated 3BD cases have become even more rare; 
phasing general n-beam cases (n > 3) is not feasible at 
the moment due to theoretical deficiencies. Therefore, it 
is important to mention that, regarding complex crystals, 
the usefulness of PPC is quite limited when compared to 
the available phasing procedures. Nevertheless, there are 
researches focused on developing and optimizing exper- 
imental data collection procedures for PPG^. On the 
other hand, the 3BD experiments offer an unique op- 
portunity for accurate determination of triplet phase in- 
variants, and consequently, for studying crystalline struc- 
tures via measurements of these invariants. For example, 
depending on the achieved experimental accuracy, elec- 
tron density of chemical bonding charges^ or even dis- 
tortion of molecules under applied electric field can be 
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FIG. 1: Experimental (open circles) and simulated (solid 
lines) (/5-scans of the 260/II2/I52 three-beam diffraction in 
a KDP crystal taken at different polarization angle, x (right- 
hand side of each scan). [001] is the reference direction (</p = 0, 
see inset). X-ray photon energy is 7482eV, and further exper- 
imental details can be found elsewhersSi. The intensity scale 
is linear, but for visualization purposes the ordinates of some 
scans are shifted from their actual values, given at the left 
(in cps). The ip-scan at x = 16° (gray scan) was mistakenly 
collected at the shoulder of the 260 reflection {Auj = 0.003°, 
30% of the FWHM = 0.01°). The flexibility of the fitting 
equation, Eq. Q, to reproduce these (p-scans is exploited in 
Fig.Hb). 



investigated by monitoring a few triplet phases. Note 
that each triplet phase is an absolute value since it al- 
ready is the phase differences between two diffracted X- 
ray waves, and not a relative quantity such as obtained 
in peak position or intensity measurements. 

This work has been motivated by our desired of devel- 
oping at LNLS a systematic and practical procedure for 
determining phase invariants with good accuracy. Exper- 
imental and analytical procedures are still to be improved 
to push phase measurements from the state-of-art to rou- 
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tinely and automatic phasing procedures; otherwise it 
will be very difficult to non-expert users to take advan- 
tages of the new possibilities offered by measuring this 
physical quantity. Data collection procedures are already 
proposed^, and undergoing improvement^, but the ac- 
tual challenger is the data analysis procedure^. Here, we 
outline the underneath principals for developing an au- 
tomatic procedure to extract accurate phase values from 
3BD interference profiles. A general systematic proce- 
dure is demonstrated, in practice, by analyzing 3BD in- 
tensity data from a KDP crystal, and the major sources 



where -Da — Dqva and Dbc{^) — D[jRf{Lp)e^^'^ v-qg 
are the amplitudes of the primary and secondary elec- 
tric displacement wavefields generated by the primary 
reflection. A, and by the detour reflection (also known as 
Umweg reflection) formed by two consecutive reflections, 
B and C. R stands for maximum amplitude ratio of these 
waves, and v^c are polarization factors for linearly 
polarized incident radiation. St is the phase difference 
between these waves, which is the triplet phase invari- 
ant. A gaussian convolution, G{u) with FWHM — wq 
and Mo = ±2.5wG, is necessary to account for the in- 
strumental width WG- f ((ys) = ws/['2{(p — ^o) — iivg] is a 
line proflle function (FWHM — w, ws — ±w) describing 
the intrinsic 3BD proflle as a function of the azimuthal 
rotation angle ip. b and ^ are related to energy balance 
mechanisms among the diffracted beams and crystalline 
imperfections, respectively^. 

Essentially, the analytical problem in accurate phase 
determination resides on how to adjust the vector of pa- 
rameters, p = [w, R,^,b,ipQ,WG], without compromising 
the extracted values for 6t- Here, a simple and fast evo- 
lutionary algorithm (DEA)» has been used for fltting the 
experimental profiles where the improvements of the fit- 
tings are guided by a mean-absolute deviation function, 
E{p). The basic strategy is then to find out the mini- 
mum of E{p) as a function of St, i.e. Eq{St), while p 
is kept within reasonable ranges of allowed values. The 
minimum of the Eq{St) curve, dEa/dST — 0, provide the 
experimental value for 6t- 



III. RESULTS AND DISCUSSIONS 

Fig-fflshows set of 3BD data collected at Brazilian Syn- 
chrotron Light Laboratory (LNLS) with the polarimeter- 
like diffractometer described elsewhere^. It is composed 
of several </3-scans taken at different polarization angles 
X, as indicated in Fig. 



of errors are pointed out. 



II. THEORETICAL BASIS 

In general, the 3BD intensity profiles are dominated 
by the interference of two diffracted waves. It leads to 
a relatively simple parametric equation that can be used 
to fit most of the experimental intensity profiles and to 
extract the phase values. It is given by^ 



(1) 

I 

Instrumental broadening effects on the interference 
profiles, as illustratively shown in Fig. [2Ia), can reduce 
accuracy when combined with the uncertainty of the R 
parameter, which is in fact the major source of inaccu- 
racy, as demonstrated in Fig.|2Ib). The £'o((5t) curves in 
Fig.|2Jb) is just showing that it is not possible to extract 
an accurate value of £^o('5t) from a single ip-scan when R 
is unknown. 

The best strategy, that we could elaborate, for ac- 
curate determination of triplet phases is composing 
polarization-dependent sets of azimuthal scans, as the 
one in Fig. ^ and then, search for the value of R that 
provides dEo/dSr = as close as possible of a same 
(5t value. Here this search strategy has been applied in 
two sets of azimuthal scans: a simulated one that is free 
of instrumental effects such as statistic noise and sample 
misalignments, and another that is the experimental data 
in Fig.Q] The Eo{6t) curves of the simulated (/?-scans for 
several values of R are shown in Fig.lSJa) while Fig.lSJb) 
shows the respective Eo{St) curves for the experimental 
data. 



IV. CONCLUSIONS 

The data analyses presented here have demonstrated 
that systematic and reliable phasing procedures are fea- 
sible. However, accuracy can be improved by optimizing 
the incident X-ray beam optics regarding energy resolu- 
tion and angular divergences, mainly in the horizontal 
plane. A good instrumental precision is also required as 
well as low noise in the intensity data. 
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FIG. 2; (a) Simulated instrumental broadening effects on (/j-scans. Simulation parameters used into Eq. (0: St = —2.6°, 
X = 32°, and p = [0.0012°, 1.0, 0.8, 0.0, 67.683°, wq] where the instrumental width values, wq, are indicated by arrows. 
A</5 = — ipo. (b) Theoretical accuracy in phase measurements as a function of the instrumental width iog, and amplitude 
ratio R. The Eo[5t) curves were obtained by fitting the profiles in (a) with wg ~ 0.001° (open circles) and wg ~ 0.006° 
(closed circles). The fittings have been carried out by the DEA within the allowed ranges: p = [0.0008° : 0.0012°, R, 0.2 : 
1.0, 0.0, 1^0 ±0.012°, 0.001° : 0.007°] where R = 1.0 (black lines) or 71 = [0.6 : 1.4] (gray lines). Definition on the dEo/dS-r = 
position gives the accuracy on St. 




FIG. 3: Absolute-mean deviation as a function of 5t, Eo{5t), obtained for (a) the simulated scans and (b) the experimental 
scans in Fig. Q All curves are normalized by its minimum value and add to an integer for better visualization. The curves with 
minima equal to 1, 2, 3, 4, 5, and 6 correspond to those scans with x = 8°, 12°, 16°, 20°, 24°, and 32°, respectively. Allowed 
range is p = [0.0010° : 0.0014°, R, 0.0 : 1.0, 0.0 : 3^^^ (fio ± 0.012°, 0.001° : 0.006°] where the R values or ranges are shown in 
the figure for each case, and changes the upper limit of the b range with the polarization angle; here = sin'^ x- 
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